Fractional spin textures in the frustrated magnet SrCr 9p Gai2 9p O 
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We consider the archetypal frustrated antiferromagnet SrCr9 P Gai2-g p Oi9 in its well-known spin- 
liquid state, and demonstrate that a Cr 3+ spin 5* = 3/2 ion in direct proximity to a pair of vacancies 
(in disordered p < 1 samples) is cloaked by a spatially extended spin texture that encodes the 
correlations of the parent spin-liquid. In this spin-liquid regime, the combined object has a magnetic 
response identical to a classical spin of length S/2 — 3/4, which dominates over the small intrinsic 
susceptibility of the pure system. This fractional-spin texture leaves an unmistakable imprint on the 
measured 71 Ga nuclear magnetic resonance (NMR) lineshapes, which we compute using Monte-Carlo 
simulations and compare with experimental data. 
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SrCrg p Gai2-9pOi9 (SCGO) is a remarkable magnetic 
material which does not show any signs of magnetic 
ordering even at extremely low temperatures T ~ 
©cvk/100, where @cw ~ 500 K is the so-called Curie- 
Weiss temperature at which mean-field theory predicts 
magnetic ordering of its corner-sharing network (Fig []} 
of antiferromagnetically coupled Cr 3+ S — 3/2 moments. 
The original observationyj-yj of this broad regime of spin- 
liquid behaviour in SCGO led to an intensive study of 
a series of samples with varying density x = 1 — p of 
vacancies in the Cr 3+ spin network, including NMR|4], 
susceptibility |5(, and //SR[6] experiments. Thus, SCGO 
is among the best known and most systematically studied 
candidates for spin-liquid behaviour in frustrated mag- 
nets. 

One key observation, due to Schiffer and Daruka, 
was the presence of a paramagnetic "Curie tail" in the 
low temperature uniform susceptibility of these com- 
pounds, well modeled by a defect contribution Xdofcct = 
Cd/T that dominates over the intrinsic susceptibility 
Xintrinsic ~ C\ / (T + Q C w) for T < QcW, the Curie con- 
stant Cd in this "two-population" phenomenology was 
associated with a population of paramagnetic objects 
dubbed "orphan spins" [5|. Another important obser- 
vation, due originally to Limot et. al.4], is the appar- 
ently symmetric NMR line broadening AH that scales 
as AH oc x/T for not-too-low x and T; this was in 
turn interpreted phenomenologically as a signature of a 
disorder-induced spin-texture — a short-ranged oscillating 
spin density profile induced by some lattice defects. 

On the theoretical side, Berlinksy and one of us used 
a "single-unit approximation" \T\ to predict that "defec- 
tive" simplices (corner-sharing units), in which all but 
one spin has been substituted for with non-magnetic im- 
purities, must give rise to Curie tails in the low tem- 
perature susceptibility of isotropic classical spin-5 1 anti- 
ferromagnets on corner-sharing lattices such as SCGO. 
In related work, Henley used classical energy minimiza- 
tion considerations at T — to argue that an infinitesi- 



mal magnetic field applied to such a system with a sin- 
gle defective simplex should induce a spin-texture with 
saturation magnetization S , /2[8j. This suggested that 
the phenomenological orphan spins [5j and oscillating spin 
textures [4| invoked earlier are related to the presence of 
such defective simplices, with the lone spins on such de- 
fective simplices providing a microscopic basis for the 
phenomenological orphan spin population of Ref S Al- 
though this orphan-texture complex (comprising the or- 
phan spin on a defective simplex and its surrounding spin 
texture) has thus been implicated in some of the most 
intriguing experimental observations on SCGO, a funda- 
mental understanding of it has been lacking. 

Here, we develop a quantitatively accurate analyti- 
cal theory that provides a full characterization of this 
orphan-texture complex by accounting for both entropic 
and energetic effects in the spin-liquid regime of low tem- 
peratures (T = kfjT/JS 2 <C 1) and low magnetic field h 
(h = gLfiFjh/JS <C 1), but arbitrary h/T. In this regime, 
we find that the orphan spin magnetisation is equal to the 
magnetisation of a spin S in a field h/2. This is accom- 
panied by an extended spin-texture, which is scale-free 
at T — but acquires a finite extent as T is increased. 
We determine the intricate pattern of spin correlations 
in the texture and demonstrate that its net magnetisa- 
tion cancels off half of the orphan spin's moment, thus 
giving rise to an orphan-texture complex that behaves as 
a classical spin S/2 in field h, with a susceptibility given 
by 

Xs/2(T) = (9L^ B S/2) 2 /ik B T 

This provides a particularly dramatic instance of frac- 
tionalisation of spins in a simple classical system. 

In remarkable correspondence with experiment [J], we 
find that these extended, fractional-spin orphan-texture 
complexes show up prominently in SCGO at not-too-low 
x = 1 — p as a large, nearly symmetric low tempera- 
ture broadening AH oc x/T of our predictions for the 
n Ga(4f) NMR lineQ associated with Ga nuclei at the 
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FIG. 1: a) The Cr atoms in SCGO form a lattice made up of 
kagome bilayers separated from each other by a layer of iso- 
lated dimers consisting of pairs of Cr atoms. Each kagome bi- 
layer is a corner sharing arrangement of tetrahedra and trian- 
gles made up of two Kagome lattices that are coupled to each 
other through "apical" Cr sites shared between up-pointing 
and down-pointing tetrahedra. Links between near-neighbour 
Cr sites in each bilayer represent a Heisenberg exchange cou- 
pling J = 80 K between neighbouring Cr 3+ spins, while links 
in the isolated dimer layer represent a Heisenberg exchange 
coupling J' = 216 K between the two Cr 3+ ions that con- 
stitute each pair. Two vacancies in a triangle (green circles) 
leave behind an orphan spin, b) The 12 Cr 3+ sites (black) 
that are hyperfine coupled to a given Ga(4f) site (red) in the 
SCGO lattice. 



so-called (4f) crystallographic position in the SCGO lat- 
tice (Fig [lj . We therefore focus below on this particular 
case, although our analytical low temperature, low field 
results apply more generally to orphan-texture complexes 
surrounding isolated defective simplices in kagome, py- 
rochlore and SCGO lattices. 

SCGO is described by the classical Hamiltonian: 
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where Si are classical length-,!? spins (S = 3/2 for SCGO) 
with saturation moment qlS in units of the Bohr mag- 
neton u B (<?£ = 2 for SCGO), J is the nearest- neighbour 
Heisenberg exchange coupling (J ~ 80K for SCGO 4]), 
and h is the external field. El refers to the tetrahedra 
that consisting of a triangle in a kagome layer attached 
to an apical spin in the triangular layer (Fig [I}, while A 
refers to those triangles in the kagome layer which are 
not associated with an apical spin. 

Here, we model the pure p — 1 system in terms of 
an effective free energy functional T that incorporates 
the energetics of the antifcrromagnetic interactions J on 
an equal footing with entropic effects of thermal fluctua- 
tions: 
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The unconstrained effective field fa serves as a surro- 
gate for the microscopic fixed- length spin variables Si, 
with the statistical weight of a field configuration {fa} 



being proportional to exp(— J-/T). H({fa}) is the clas- 
sical Hamiltonian of the system now written in terms of 
<f>i, and the phenomenological stiffness constants pi are 
fixed by requiring that the mean length of fa equal S 
within the effective theory: (fa 2 )? = S 2 \$. For a back- 
ground to this approach for pure systems, see the review 
by Henley 10]. 

To incorporate disorder effects in diluted p < 1 sam- 
ples, we assume that the stiffness constants do not change 
significantly from their pure values, but extend the ef- 
fective theory in two important ways: First, we model 
vacancies in the lattice by setting fa to zero on vacancy 
sites. Secondly and more crucially, as the fixed length na- 
ture of an orphan spin on a defective simplex is expected 
to play a central role, we retain it as a microscopic length- 
s' spin and do not introduce an effective field variable at 
such sites. 

We now use this effective theory to analyze the SCGO 
magnet with a single defective triangle, i.e. two vacancies 
on one triangle of the SCGO lattice (Fig[T}. In this case, 
the effective theory reduces to a length-S orphan spin Sn 
(n 2 = 1) coupled to a constrained Gaussian theory for 
fa (with constraints fa = at the two vacancy sites). 
Focusing first on the orphan spin Sn in this defective 
triangle by integrating out the <j> fields, we find[9( that the 
exchange field from the surrounding spin liquid "screens" 
exactly half the external magnetic field on the orphan 
spin in the low temperature, low field spin liquid regime, 
yielding a spin S variable that "sees" a magnetic field 
h/2, and therefore develops a polarization equal to that 
of a free classical spin S in a field h/2. This striking 
prediction is fully confirmed by Monte-Carlo studies[9| 
of the classical SCGO magnet with one defective triangle 
(see Fig [5] (a)), which also reveal that this prediction is 
surprisingly robust, remaining accurate for temperatures 
as high asT~ O.IJS 2 for this example. 

Next, we turn to a detailed description of the extended 
spin texture that cloaks this orphan spin in the spin- 
liquid regime. In addition to the uniform external field 
h = hz that acts on all the (f>i in the constrained Gaussian 
action, this orphan spin polarization also gives rise to a 
local exchange field that acts in the z direction on fa 
at the two undiluted sites adjacent to the orphan spin. 
The extended spin texture surrounding the orphan spin 
is modeled within the effective theory by calculating the 
response (<j>* )f to these fields[9|]. At T = 0, the computed 
texture's envelope decays as l/\f\ away from the orphan 
spin, while the overall scale is set by the orphan spin's 
saturation magnetization. At finite-T and small fields in 
the spin-liquid regime, the power law envelope is cut off 
by a thermally introduced finite correlation length £ ~ 
1/y/T, endowing it with an effective size £ d ~ 1/T d / 2 in d 
dimensions; in our d — 2 example, this gives a spatial size 
~ 1/T scaling identically with the overall magnetization 
scale of the texture, which is set by the orphan spin's 
susceptibility ~ 1/T. These predictions are compared 
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FIG. 2: a) The orphan spin in field h develops magnetization (symbols) that matches that of an isolated spin S in field h/2 
(full line) in the spin liquid regime. Also shown (dotted curve) is the magnetization expected for an isolated spin 5* in field h. 
b) The orphan induced spin texture is shown along two cuts on the lattice. Numerical data is shown as points with statistical 
error bars and the effective theory result as lines in the main plot. The spin texture shown in red (blue) in main plot is along 
the cut shown in red (blue) in the inset. The orphan spin is denoted by an arrow in the inset, c) The impurity magnetization 
(defined as the difference in magnetic response between the diluted and the pure system) of a single orphan-texture complex 
(symbols) compared with Bs/ 2 (h,T) (solid curve), the asymptotic prediction deep in the spin-liquid regime. Bs/ 2 {h,T) is the 
magnetization curve of a classical spin 5/2. Note the much faster convergence of numerical data to asymptotic predictions in 
the linear (small hS/T) part of the magnetization curve as compared with the non- linear regime at large hS/T. Inset: The 
corresponding impurity susceptibility (defined in text) matches the susceptibility of a free classical spin 5/2 up to remarkably 
high temperatures T ~ 0.IJ5 2 . All Monte-Carlo data exhibited in a) was extrapolated to the thermodynamic limit using a 
sequence of sizes up to L — 36 (with 7L 2 sites), while that in c) was extrapolated using a sequence of sizes up to L — 20. 



with Monte-Carlo simulation[9| results in Fig [5] (b) for 
the case of a single defective triangle in SCGO, and the 
agreement is remarkably good. 

Although the computed texture decays slowly (as l/|r| 
at T = 0), its oscillations conspire to reduce its net 
moment so that it cancels out precisely half of the or- 
phan spin's moment, endowing the orphan-texture com- 
plex with a net "impurity magnetization" (excess mag- 
netization over and above that of a pure system in the 
same external field) equal to the magnetization of a clas- 
sical spin S/2 in field h. The corresponding "impurity 
susceptibility" Ximp arising from a single orphan-texture 
complex is (gL^BS/2) 2 /3fcsT — a Curie tail dominating 
the low-T response. For the case of a single defective 
triangle in SCGO, our Monte Carlo simulations validate 
this striking prediction up to surprisingly high tempera- 
tures of order T ~ OAJS 2 (Fig [2] (c)). 

Fortunately, this novel physics leaves its imprint on the 
Ga(4f ) NMR line, which is a sensitive probe of the distri- 
bution of the twelve Cr spin polarizations that are hyper- 
fine coupled to each Ga(4f) nucleus in SCGO (FigQ](b)). 
We calculate this distribution in diluted SCGO lattices 
with uncorrelated site dilution probability x — 1 — p[9(, 
using the experimentally known values [4j of these hyper- 
fine couplings, and Monte-Carlo simulations [9j to fully 
account for the physics of a non-zero density of orphan 
spins in order to obtain quantitatively accurate results. 

As is clear from the example at x = 0.2 shown in Fig [3] 
(a) , this calculation predicts a line that is broad and ap- 
pears symmetric, reflecting the fact that the spin textures 
are staggered and involve a very large number of spins, 
making it difficult to discern the 0(1) net moment of each 



texture that endows the line with a slight bias towards 
lower magnetic fields. 

In Fig [3] (c) , we show the temperature and x depen- 
dence of the width of our predicted lines for not-too-small 
values of x and T, of greatest relevance to experiments. 
As the lineshape is not well-approximated by a Gaus- 
sian, we do not fit the line to a Gaussian, but rather use 
the definition AH = 2^/2\n(2) ((H 2 ) - (H) 2 ) 1/2 which 
reduces to the standard value for a Gaussian line but 
provides an unbiased measure of the width in more gen- 
eral cases. As is clear from Fig [3] (b), the theoreti- 
cal preditions for AH have the expected "Curie tail" 
AH ~ A(x)/T at low temperature. For uncorrelated 
site disorder, the coefficient A is expected to scale as x 2 
for asymptotically small x; however, for not-too-small x 
relevant to experiment, we find that A{x) can be fit well 
by an approximately linear x dependence A(x) ~ x. 

Thus, our theory reproduces the broad, apparently 
symmetric lines seen in experiments, with broadening 
scaling as AH ~ x/T for not-too-low x and T. Although 
these experimental facts are fully reproduced, we caution 
that our results do not provide a fully quantitative ex- 
planation of the NMR data: As is clear from Fig [3] (b), 
our minimal disorder model consisting of independent va- 
cancies of the nominal experimental density significantly 
underestimates the absolute scale of the linewidth AH. 
Since isolated vacancies lead to no oscillating spin tex- 
tures scaling as 1/T|9|, this discrepancy clearly demon- 
strates that the density of orphan spins in the experi- 
mental samples is significantly higher than that expected 
from uncorrelated Ga substitution of the Cr lattice. The 
most frugal resolution is perhaps that the Ga substitu- 
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FIG. 3: Predictions for Ga(4f) NMR lines, a) 71 Ga(4f) "field-scan" NMR lines predicted within the minimal disorder model 
of uncorrelated dilution at vacancy density x — 0.2. (b) Experimental results of Ref [J] for field-scan 71 Ga(4f) NMR lines 
(dots) at dilution x = 0.19, plotted using datafiles provided by P. Mendels, and compared to theoretical predictions of the 
minimal disorder model (with y axis rescaled and x axis offset appropriately for easy comparison of linewidths and shapes of 
all curves) at vacancy density x = 0.2 (dashed curves) and x — 0.3 (full curves). Only the low-field side of the experimental 
line is compared with theory, as the high-field side of the experimental line is known to be contaminated by the presence of a 
satellite peak arising from non-stoichiometric (impurity) 71 Ga nuclei located at three different Cr sites in the SCGO lattice[J]. 
(c) Temperature and x dependence of the width AH of the theoretically predicted field-scan 71 Ga(4f) NMR lines, normalized 
by a reference field H re f = 3.12T. Solid curves in the main plot are fits of AH to A(x)/T over the experimentally relevant 
temperature range, and inset shows the approximately linear x dependence of the best-fit values of A(x) over the experimentally 
relevant range of x. The Monte-Carlo simulations were performed using a sequence of sizes ranging from L = 16 to L = 50 
(with 7L 2 sites) to eliminate finite-size effects, with on average 20 disorder realizations at each size. 



tion in the Cr kagome layer may be correlated. Addi- 
tionally, and perhaps more realistically, other sources of 
disorder, for instance random strain-induced bond ran- 
domness, could also affect the effective density of 'or- 
phan spins'. Unfortunately, it does not seem possible to 
use measurements of the impurity Curie constant Cd[5j 
to shed light on this, as uniform susceptibility measure- 
ments also pick up, apart from the signal due to Kagome 
layer orphan spin textures, a much larger "background" 
contribution from free Cr spins created when Ga impuri- 
ties substitute for Cr in the isolated dimer layer (Fig [T]) 
of the SCGO lattice. 

In conclusion, we note that these disorder induced 
orphan-texture complexes in SCGO embody several cen- 
tral themes of modern condensed matter physics: The 
emergence of new types of extended degrees of freedom, 
and their rich physics at finite density, is a common 
strand that runs through diverse examples llj such as 
Skyrmion lattices in spinor condensates|12J|. itinerant |13| 
and quantum Hall magnets 141 Il5ll , and spin textures 
in topological phases of matter [16j. Moreover, the fact 
that our textures give rise to a relatively strong magnetic 
response — their staggered l/r d ~ l magnetisation profile is 
parametrically stronger than the intrinsic l/r d spin corre- 
lations of the parent spin-liquid — shows that, like in the 
case of defects evidencing the d-wave nature of the order 
parameters in the cuprates, imperfections are among the 
best and most direct probes of exotic correlated states of 
matter. 
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